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Abstract 

We have observed 23.2±6.01q g purely-leptonic decays of Df — > /i + f /t from a sample 
of muonic one prong decay events detected in the emulsion target of Fermilab exper- 
iment E653. Using the — > 0/x + z^ yield measured previously in this experiment, 
we obtain B(D+ -> ^ + v^)/B(D+ -» = 0.16 ± 0.06 ± 0.03. In addition, we 

extract the decay constant f Dg = 194 ± 35 ± 20 ± 14 MeV. 
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Leptonic decay of charged pseudoscalar mesons is described by the annihilation of con- 
stituent quark and antiquark into a virtual W boson. The rate for this process is proportional 
to the square of the pseudoscalar decay constant /, which measures the overlap of quark 
and antiquark at zero separation. The constant / gives absolute normalizations of numerous 
heavy-flavor transitions, including mixing, semi-leptonic and non-leptonic decays. Hence it 
is highly desirable to determine the pseudoscalar decay constant more precisely in charm and 
B physics. A measurement of a purely-leptonic decay branching ratio, unlike semi-leptonic 
or non-leptonic decays, is the most reliable way to extract the pseudoscalar decay constant, 
because it does not involve any QCD corrections. 

There are three purely-leptonic decay modes for the (charge conjugate modes are 
implied throughout this paper): electronic, muonic and tauonic. Among them, the muonic 
decay is the most accessible mode to determine the decay constant fo 3 ■ The electronic decay 
is helicity suppressed by five orders of magnitude with respect to the muonic one. While 
the tauonic decay branching ratio is expected to be an order of magnitude larger than the 
muonic decay, the Q value for — > t + v t is small and hence reconstruction of this mode is 
extremely difficult because of the small decay angle in the lab frame. 

The first measurement of the pseudoscalar decay constant of the charmed strange meson, 
fo s , was carried out by the WA75 group in 1993 [0], using a transverse momentum spectrum 
of muons from leptonic decays observed in an emulsion target. They obtained a value of 
f Ds = (232 ±45 ±20 ±48) MeV. CLEO p( and BES P] also reported recent measurements 
of fj} s , and ARGUS pi] obtained fo s by model dependent calculation in 1992. All of these 
results, however, suffer from large statistical and systematic errors, and hence the present 
measured value of fo 3 is still not determined well enough to discriminate between different 
theoretical models. In this paper, we report an independent measurement of fo s from an 
analysis of the data from Fermilab experiment E653. 

The Fermilab emulsion-hybrid experiment E653 was designed to study production and 
decay of heavy flavor particles by the direct observation of decay vertices in the emulsion. 
Data for this analysis were taken in a 600GeV/c tt~ beam during the second run of Fermilab 
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E653. The trigger required both a beam particle to interact in the target and a muon 
to penetrate 3,900 g/cm 2 of absorber. This muon trigger provided an enriched sample of 
semi-muonic decays of heavy flavor particles. A total of 8.2 x 10 6 events, corresponding to 
2.5 x 10 8 interactions, were recorded during this run. Since the details of this experiment 
were described in another paper P] , only a brief description of each detector element is given 
here. 

The detector consisted of the emulsion target, tracking detector, spectrometer, and muon 
detection system. Because of its unsurpassed spatial resolution, nuclear emulsion was chosen 
as an active target, and it was used both as a primary interaction target and as a decay vertex 
detector. The electronic detectors downstream of the target select events for the emulsion 
scanning and predict the interaction vertices and each track position in the emulsion. Tracks 
from interactions and decays in the emulsion were measured in the vertex silicon strip 
detectors (VSSDs). The VSSDs consisted of 18 planes arranged in 6 groups of three readout 
coordinates, and the center part of the detector achieved a resolution transverse to the 
beam of 8.8/xm. This detector also served as the upstream arm of the spectrometer. The 
spectrometer gave a momentum resolution of cr p /p = <J (0.01) 2 + (0.00023p) 2 (p in GeV/c). 
The most down stream element of the E653 detector array was the muon identification 
system, which remeasured the direction and momenta of muon candidates after the absorber 
with drift chambers and an iron toroid. 

To enrich the relative yield of events with heavy flavor particle decays and to reduce 
the emulsion scanning load, the transverse momentum of the triggered muon with respect 
to the beam axis, Pt^ Beam, was required to be greater than O.SGeV/c. In total, 94,000 
events passed this criterion from a initial sample of 8.2 x 10 6 events. This reduced sample 
contained a significant number of background events from feedthrough muons. Pions or 
kaons could be misidentified as muons because they decayed muonically in flight with small 
angle deflections or from punch through in the iron absorber. These events were rejected 
using the method described below. 

In the emulsion analysis, the primary interaction vertices were visually located in the 
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emulsion by a semi-automatic scanning microscope system using position predictions from 
the VSSDs. The efficiency for finding primary vertices was determined to be greater than 
99%. If the angle of an emulsion track from the primary vertex matched to within 1 mrad 
of the reconstructed muon track, the muon was identified as a feedthrough muon and the 
event was rejected. For the remaining events, a more detailed emulsion scanning technique 
was employed to find the origin of the muon. 

Different decay search methods were applied for charged and neutral decay vertices. 
To find charged decays, all emulsion tracks observed at the primary vertex which did not 
match to any track reconstructed by the VSSDs were followed towards the most downstream 
emulsion plates. To find neutral decays, the tracks reconstructed by the VSSDs with no 
corresponding emulsion track observed at the primary vertex, were located at the most 
downstream plate of the emulsion target. Once found, these tracks were followed upstream 
through the emulsion stack to their origin within the emulsion target to detect neutral 
decays. The detailed event location and vertex search methods are described elsewhere P]. 

As a result of the above analysis, 1,193 charmed particle candidates were identified: 565 
KINK (charged one prong) decays, 536 VEE (neutral two prong) decays, and 92 three or four 
prong decays. Further cuts on the decay length (> 500//m) and the transverse momentum 
of the muon with respect to the parent particle direction (-Pz> > 0.28GeV/c) were applied 
to the KINK and VEE decays to improve position resolution measurement and to reject 
strange particle decays. Only VEEs with their vertices found in the emulsion were retained 
to avoid the large Pt^ error due to the intrinsically larger position errors of the VSSDs. The 
offline selection criteria and cuts are summarized in Table |l|. 

TABLE 

Figure [I] shows the Pt^ distributions for 531 KINKs and 276 VEEs which survived the 
above selection criteria. An excess of 19 events can be seen in the region (Pt/i > 0.8Q5GeV/c) 
exceeding the kinematic limit of the D + — > K^^v^ decay for the KINK events, while no 
such signal is seen for the VEE events. These large Pt^ events can be interpreted as the 
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fully leptonic decay Df — > /U + z/ M . 

In addition to — ► this sample includes background events. The largest source 
of background is the muon from the decay D + — > p^v^. The inclusive production cross 
section for D + is a factor of two larger than D s . However, the branching fraction for this 
decay is considerably smaller than Df — > /i + z/ M , because it is suppressed by the ratio of 
CKM matrix-element, | V cd | 2 / | V cs \ 2 . A muon from the cascade decay, — > t + u t ; 
t + — > h+v^Ut is another source of background. Since the decay angle of the first decay point 
is extremely small (^1-2 mrad), it is often not detectable by our ordinary emulsion scanning 
method [i]. Therefore, only the second decay r + — > /i + z/ M z7 T is observable and therefore will 
appear in the emulsion as a muonic KINK. A final possible background is a feedthrough 
muon from the daughter track of a hadronic charm decay. 

We used a Monte Carlo simulation to estimate the total number of background events. 
Charged and neutral charmed particles were generated, using the production distribution 
parametrization (fa/dxpdP^ oc (1— | xp \) n • exp(-bP^) with n = 4.25 ± 0.24 ± 0.23 and 
b = 0.76 ± 0.03 ± 0.03(GeV/c)~ 2 P]. The relative numbers of D + ,D° and Df mesons were 
generated according to the the cross section values measured in this experiment lijfll. The 
baryon was generated using the cross section value obtained from other experiments [0]. 
The cross section numbers employed in the Monte Carlo are summarized in Table |^. Decays 
of charmed particles were simulated by JETSET 7.4 with known branching fractions [0], 
except that the ratio /d+//d s = 0.90 [0] was used for B(D + — > /U + z/^). 

TABLE | 

The number of — > was then extracted from a fit to the Pt^ distribution of the 
muonic KINK sample. In this fit, the background normalizations and shapes were fixed by 
the Monte Carlo, while the normalization for the Df — > yU + z^ M decay was allowed to float. 
Fig. [I] a) shows the fit result for the muonic KINK decays. The points with error bars are 
the experimental data, while the solid, dotted and dashed histograms represent the final fit 
result, the sum of the backgrounds and the — > signal, respectively. The fit is in 
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a very good agreement with the data, with / D.F. is 1.029. The fit yielded 23.2 ± 6.0lo;g 
D+ — > /U+z/^ decays. The first error is the statistical and the second one is the systematic, 
involving uncertainty of production distribution. The detection efficiency for these decays is 
14.6%, determined mainly from a Monte Carlo calculation. It includes efficiency factors for 
trigger and reconstruction (0.817), offline selection cut (0.416), fraction of sample scanned 
(0.570), fiducial and location losses (0.956), final selection cuts (0.855), and scanning effi- 
ciency in emulsion (0.925). The overall efficiency is large and well- understood because of 
the simple kink topology, the weak event-selection cuts, and the high scanning efficiency. 
The number of the background events from D + — > /i + z/ M and — > t + u t ; t + — > ^i + v^i7 T in 
the high Pt^ region were estimated to be 4.5 and 0.28 respectively. The final background, 
feedthrough muons from hadronic decays, was determined to be negligible. As a systematic 
check, the Monte Carlo prediction for the VEE decays was compared with the data as shown 
in Fig. |l]b). They are in a good agreement, which gives further confidence in our analysis. 

FIGURE 

In a previous publication P], we reported on 18.7±4.9lg;7 D+ — > events measured 

in E653 without emulsion information with an efficiency of 1.86%. Combining these results, 
we can calculate the ratio B(D+ — ► n + v^) / B(D+ — > 0/i+z/^) using the following equation: 



0.16 ±0.06 ±0.03, (1) 



where e and N are the detection efficiency and the extracted yield for each decay channel 
respectively. The first error is the statistical and the second one is the systematic. Using 
the measured value of B(D+ -> <f>l+i^) = 1.88 ± 0.29% [0], we obtain, 



B(Df -> = (0.30 ± 0.12 ± 0.06 ± 0.05)% (2) 

The first two errors are the statistical and systematic errors of our analysis, involving un- 
certainty of production cross section and its distribution, and the third one reflects the 
uncertainty in the branching fraction for — > 0/ + z/^. Note that this normalization does 



not require knowledge of the cross section. Furthermore, it is quite insensitive to the 
production parameters n and b because the muon Pj- relative to the beam direction in 
— ► /U + z/ M is dominated by the decay Pt of the muon rather than by production Pt of the 

The pseudoscalar decay constant fn a can be extracted from the following equation, 
B{Dt - ^) = ^-f D T D m D ml \ V cs | 2 (1 - (3) 

<57T " l D s 

where Gp is the Fermi constant, Tp> g and mp> g are the mean lifetime and mass of the D+, 
m M is the mass of muon, and V cs is the Kobayashi-Maskawa matrix element. Using Tp> s = 
4.67 ± 0.17 x 10~ 13 s, m Ds = 1968.5 ± 0.6 MeV/c 2 , | V cs |= 0.9745 ± 0.0007 P], and our 
result for B(Df — > we obtain, 

f Da = 194 ± 35 ± 20 ± 14 MeV (4) 

The first two errors are the statistical and systematic errors of our analysis, and the third 
one reflects the uncertainty in the branching fraction for Df — > <pl + vi. 

Table |3| compares the E653 result for fr> s with those from other experiments. Agreement 
with WA75 pt| is excellent. Our result is also consistent at the ~10% confidence level with 
the lager values of fo s from CLEO [!] and BES P]. Theoretical predictions span 
the range 130 to 350 MeV . Recent calculations using QCD sum rule pi], independent quark 
model |@], and lattice P,0| techniques are in good agreement with our result. 

TABLE | 

In conclusion, E653 has measured the ratio of the muonic branching ratios for the D s ; 
B(D+ -> fi + v^)/B(Df -> <M + ^) = 0.16 ± 0.06 ± 0.03. With the latest value for B(D+ -> 
# + ^), we obtain B{D+ -> ^v^) = (0.30 ± 0.12 ± 0.06 ± 0.05)%. This yields a value of 
f Da = 194 ± 35 ± 20 ± 14 MeV. 

We gratefully acknowledge the efforts of the Fermi National Accelerator Laboratory 
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FIGURES 

Fig. 1. Decay Pt^ distributions of a) 531 KINK (muonic one prong) events, and b) 276 VEE 
(muonic two prong) events. The points with error bars are the experimental data, and the his- 
tograms are the Monte Carlo simulations. The vertical dotted lines show the decay Pt^ limits of 
the main semi-leptonic modes. In (a), the solid, dotted and dashed histograms represent the final 
fit result, the sum of the backgrounds, and the Df — > /U + f M signal, respectively. In (b), the solid 
line represents the inclusive spectrum of muons from D° decay. 
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TABLES 



Table 1. Summary of offline selection criterion and cut requirements. 



Offline selection 
Flight length cut 
Decay Pt^ cut 

Vertex position cut (only for VEE) 



Pr^Beam > 0.8 GeV/c 

FL > 500 fim 

P T ^ > 0.28 GeV/c 

Decay vertex in emulsion bulk 



Table 2. Cross sections used for Monte Carlo background calculations. 



a(D ± ;x F > 0) 8.66 ± 0.46 ± 1.96 fib/nucleon | 

a(D°; x F > 0) 22.05 ± 1.37 ± 4.82 fib/nucleon 

a(Df; x F > 0) 5.1 ± 1.3 ± 1.3 fib/nucleon 

a(Af; x F > 0) 15.3 ± 1.7 ± 2.6 fib/nucleon f 



Table 3. Summary of published measurement of pseudoscalar decay constant fo s ■ 



Experiment fo s (MeV) 

E653 (This work) 194 ± 35 ± 20 ± 14 

WA75 [§] (1993) 238 ± 47 ± 21 ± 43 

CLEO H (1994) 344 ± 37 ± 52 ± 42 

BES [I] (1995) 4301^ ± 40 
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Abstract 

We have observed 23.2±6.0^q q purely-leptonic decays of Df — > i-i + v ll from a sample 
of muonic one prong decay events detected in the emulsion target of Fermilab exper- 
iment E653. Using the Df — > ^^v^ yield measured previously in this experiment, 
we obtain B(D+ -> n + v^)/ B(D+ -> c^/i+z^) = 0.16 ± 0.06 ± 0.03. In addition, we 
extract the decay constant f Ds = 194 ± 35 ± 20 ± 14 MeV . 
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Leptonic decay of charged pseudo-scalar mesons is described by the annihilation of con- 
stituent quark and antiquark into a virtual W boson. The rate for this process is proportional 
to the square of the pseudoscalar decay constant /, which measures the overlap of quark and 
antiquark at zero separation. The constant / gives absolute normalizations of numerous 
heavy-flavor transitions, including mixing, semi-leptonic and non-leptonic decays. Hence it 
is highly desirable to determine the pseudoscalar decay constant more precisely in charm and 
B physics. A measurement of a purely-leptonic decay branching ratio, unlike semi-leptonic 
or non-leptonic decays, is the most reliable way to extract the pseudoscalar decay constant, 
because it does not involve any QCD corrections. 

There are three purely-leptonic decay modes for the Df (charge conjugate modes are 
implied throughout this paper): electronic, muonic and tauonic. Among them, the muonic 
decay is the most accessible mode to determine the decay constant fu s . The electronic decay 
is helicity suppressed by five orders of magnitude with respect to the muonic one. While 
the tauonic decay branching ratio is expected to be an order of magnitude larger than the 
muonic decay, the Q value for Df — > t + u t is small and hence reconstruction of this mode is 
extremely difficult because of the small decay angle in the lab frame. 

The first measurement of the pseudoscalar decay constant of the charmed strange meson, 
/d s , was carried out by the WA75 group in 1993 f 1 ], using a transverse momentum spectrum 
of muons from Df leptonic decays observed in an emulsion target. They obtained a value of 
f Ds = (232 ± 45 ± 20 ± 48) MeV. CLEO [ 2 ] and BES [ 3 ] also reported recent measurements 
of fr> s} and ARGUS [ 4 ] obtained fr> s by model dependent calculation in 1992. All of these 
results, however, suffer from large statistical and systematic errors, and hence the present 
measured value of fr> s is still not determined well enough to discriminate between different 
theoretical models. In this paper, we report an independent measurement of fr> s from an 
analysis of the data from Fermilab experiment E653. 

The Fermilab emulsion-hybrid experiment E653 was designed to study production and 
decay of heavy flavor particles by the direct observation of decay vertices in the emulsion. 
Data for this analysis were taken in a 600GeV/ c ir~ beam during the second run of Fermilab 
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E653. The trigger required both a beam particle to interact in the target and a muon 
to penetrate 3,900 g / 'cm 2 of absorber. This muon trigger provided an enriched sample of 
semi-muonic decays of heavy flavor particles. A total of 8.2 X 10 6 events, corresponding to 
2.5 X 10 8 interactions, were recorded during this run. Since the details of this experiment 
were described in another paper [ 5 ], only a brief description of each detector element is given 
here. 

The detector consisted of the emulsion target, tracking detector, spectrometer, and muon 
detection system. Because of its unsurpassed spatial resolution, nuclear emulsion was chosen 
as an active target, and it was used both as a primary interaction target and as a decay 
vertex detector. The electronic detectors downstream of the target select events for the 
emulsion scanning and predict the interaction vertices and each track position in the emulsion. 
Tracks from interactions and decays in the emulsion were measured in the vertex silicon strip 
detectors (VSSDs). The VSSDs consisted of 18 planes arranged in 6 groups of three readout 
coordinates, and the center part of the detector achieved a resolution transverse to the beam of 
S.Sfim. This detector also served as the upstream arm of the spectrometer. The spectrometer 
gave a momentum resolution of (7, p /p = ^(O.Ol) 2 + (0.00023p) 2 (p in GeV/c). The most 
down stream element of the E653 detector array was the muon identification system, which 
remeasured the direction and momenta of muon candidates after the absorber with drift 
chambers and an iron toroid. 

To enrich the relative yield of events with heavy flavor particle decays and to reduce the 
emulsion scanning load, the transverse momentum of the triggered muon with respect to 
the beam axis, Py^ Beam, was required to be greater than O.SGeV/ c. In total, 94,000 events 
passed this criterion from a initial sample of 8.2 X 10 6 events. This reduced sample contained 
a significant number of background events from feedthrough muons. Pions or kaons could be 
misidentified as muons because they decayed muonically in flight with small angle deflections 
or from punch through in the iron absorber. These events were rejected using the method 
described below. 

In the emulsion analysis, the primary interaction vertices were visually located in the 
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emulsion by a semi-automatic scanning microscope system using position predictions from 
the VSSDs. The efficiency for finding primary vertices was determined to be greater than 
99%. If the angle of an emulsion track from the primary vertex matched to within 1 mrad of 
the reconstructed muon track, the muon was identified as a feedthrough muon and the event 
was rejected. For the remaining events, a more detailed emulsion scanning technique was 
employed to find the origin of the muon. 

Different decay search methods were applied for charged and neutral decay vertices. 
To find charged decays, all emulsion tracks observed at the primary vertex which did not 
match to any track reconstructed by the VSSDs were followed towards the most downstream 
emulsion plates. To find neutral decays, the tracks reconstructed by the VSSDs with no 
corresponding emulsion track observed at the primary vertex, were located at the most 
downstream plate of the emulsion target. Once found, these tracks were followed upstream 
through the emulsion stack to their origin within the emulsion target to detect neutral decays. 
The detailed event location and vertex search methods are described elsewhere [ 6 ]. 

As a result of the above analysis, 1,193 charmed particle candidates were identified: 565 
KINK (charged one prong) decays, 536 VEE (neutral two prong) decays, and 92 three or four 
prong decays. Further cuts on the decay length (> 500/im) and the transverse momentum of 
the muon with respect to the parent particle direction (Pz> > 0.28GeV/c) were applied to 
the KINK and VEE decays to improve position resolution measurement and to reject strange 
particle decays. Only VEEs with their vertices found in the emulsion were retained to avoid 
the large Py^ error due to the intrinsically larger position errors of the VSSDs. The offline 
selection criteria and cuts are summarized in Table 1. 

TABLE 1 

Figure 1 shows the Py^ distributions for 531 KINKs and 276 VEEs which survived the 
above selection criteria. An excess of 19 events can be seen in the region (Pz> > 0.865GeV r /c) 
exceeding the kinematic limit of the D + — > K°/i + u ll decay for the KINK events, while no 
such signal is seen for the VEE events. These large Py^ events can be interpreted as the 
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fully leptonic decay Df — > ji^v^. 

In addition to Df — > /i + ^u this sample includes background events. The largest source of 
background is the muon from the decay D + — > y^v^. The inclusive production cross section 
for D + is a factor of two larger than D s . However, the branching fraction for this decay is 
considerably smaller than Df — > /i + z^, because it is suppressed by the ratio of CKM matrix- 
element, | V c d | 2 / | V cs | 2 . A muon from the cascade decay, Df — > t + u t ; r + — > ji + h , ll i7 T 
is another source of background. Since the decay angle of the first decay point is extremely 
small (ssl-2 mrad), it is often not detectable by our ordinary emulsion scanning method [ 6 ]. 
Therefore, only the second decay r + — > ji + h , ll i7 T is observable and therefore will appear in the 
emulsion as a muonic KINK. A final possible background is a feedthrough muon from the 
daughter track of a hadronic charm decay. 

We used a Monte Carlo simulation to estimate the total number of background events. 
Charged and neutral charmed particles were generated, using the production distribution 
parametrization d 2 a / dxpdPj oc (1— | xp \) n ■ exp(-bPj) with n = 4.25 ± 0.24 ± 0.23 and 
b = 0.76 ± 0.03 ± 0.03(GeV r /c)~ 2 [ 8 ]. The relative numbers of D + } D° and Df mesons were 
generated according to the the cross section values measured in this experiment [ 8 , 9 , 18 ]. The 
baryon was generated using the cross section value obtained from other experiments [ 10 ]. 
The cross section numbers employed in the Monte Carlo are summarized in Table 2. Decays 
of charmed particles were simulated by JETSET 7.4 with known branching fractions [ 18 ], 
except that the ratio fo+/fD s = 0.90 [ 7 ] was used for B(D + — > 

TABLE 2 

The number of Df — > /i + ^u was then extracted from a fit to the Py^ distribution of the 
muonic KINK sample. In this fit, the background normalizations and shapes were fixed by 
the Monte Carlo, while the normalization for the Df — > /i + ^u decay was allowed to float. 
Fig. 1 a) shows the fit result for the muonic KINK decays. The points with error bars are the 
experimental data, while the solid, dotted and dashed histograms represent the final fit result, 
the sum of the backgrounds and the Df — > /i + ^u signal, respectively. The fit is in a very good 
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agreement with the data, with x 2 /D.F. is 1.029. The fit yielded 23.2 ± 6.0^:° Dj -> /i+i^ 
decays. The first error is the statistical and the second one is the systematic, involving 
uncertainty of production distribution. The detection efficiency for these decays is 14.6%, 
determined mainly from a Monte Carlo calculation. It includes efficiency factors for trigger 
and reconstruction (0.817), offline selection cut (0.416), fraction of sample scanned (0.570), 
fiducial and location losses (0.956), final selection cuts (0.855), and scanning efficiency in 
emulsion (0.925). The overall efficiency is large and well-understood because of the simple 
kink topology, the weak event-selection cuts, and the high scanning efficiency. The number 
of the background events from D + — > /i + ^u and Df — > t + u t ; t + — > ji + h , ll i7 T in the high Py^ 
region were estimated to be 4.5 and 0.28 respectively. The final background, feedthrough 
muons from hadronic decays, was determined to be negligible. As a systematic check, the 
Monte Carlo prediction for the VEE decays was compared with the data as shown in Fig. 1 
b). They are in a good agreement, which gives further confidence in our analysis. 

FIGURE 1 

In a previous publication [ 9 ], we reported on 18.7 ±4.9^o!7 — > cji^v^ events measured 
in E653 without emulsion information with an efficiency of 1.86%. Combining these results, 
we can calculate the ratio B(Df — > fi+v^)/ B(Df — > c/>/i + i/ |U ) using the following equation: 

H^I ^ f =^EL.lEL = 0.16 ± 0.06 ± 0.03, (1) 

where e and N are the detection efficiency and the extracted yield for each decay channel 
respectively. The first error is the statistical and the second one is the systematic. Using the 
measured value of B(D+ (f)l + vi) = 1.88 ± 0.29% [ 18 ], we obtain, 

B(Df n+vp) = (0.30 ± 0.12 ± 0.06 ± 0.05)% (2) 

The first two errors are the statistical and systematic errors of our analysis, involving un- 
certainty of production cross section and its distribution, and the third one reflects the 
uncertainty in the branching fraction for Df — > 4>l + fi. Note that this normalization does 
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not require knowledge of the Df cross section. Furthermore, it is quite insensitive to the 
Df production parameters n and b because the muon Py relative to the beam direction in 
Df — y /i + i/u is dominated by the decay Py of the muon rather than by production Py of the 



The pseudoscalar decay constant fr> s can be extracted from the following equation, 



where Gf is the Fermi constant, td s and my), are the mean lifetime and mass of the Df , 
is the mass of muon, and V cs is the Kobayashi-Maskawa matrix element. Using td s = 



4.67 ± 0.17 x 10" 13 s, m Ds = 1968.5 ± 0.6 MeV/c\ \ V cs | = 0.9745 ± 0.0007 [ 18 ], and our 
result for B(Df — > /i+z/^), we obtain, 



The first two errors are the statistical and systematic errors of our analysis, and the third 
one reflects the uncertainty in the branching fraction for Df — > 4>l + fi. 

Table 3 compares the E653 result for fr> s with those from other experiments. Agreement 
with WA75 f 1 ] is excellent. Our result is also consistent at the ~10% confidence level with 
the lager values of fr> s from CLEO [ 2 ] and BES [ 3 ]. Theoretical predictions [ n , 12 , 13 , 14 ] span 
the range 130 to 350 MeV . Recent calculations using QCD sum rule [ 15 ], independent quark 
model [ 16 ], and lattice [ 7 , 17 ] techniques are in good agreement with our result. 



In conclusion, E653 has measured the ratio of the muonic branching ratios for the D s ; 
B(D+ ti + u ll )/B{D+ <M + ^) = 0.16 ± 0.06 ± 0.03. With the latest value for B(D+ 
(f)fi + u^, we obtain B(D+ /i+i^) = (0.30 ± 0.12 ± 0.06 ± 0.05)%. This yields a value of 
f Ds = 194 ± 35 ± 20 ± 14 MeV. 

We gratefully acknowledge the efforts of the Fermi National Accelerator Laboratory staff 
in staging this experiment. This work was supported in part by US Department of Energy; 




(3) 



fD. 



194 ± 35 ± 20 ± 14 MeV 



(4) 



TABLE 3 
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Japan-US Cooperative Research Program for High Energy Physics; the Ministry of Educa- 
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Basic Science Research Institute Program, Ministry of Education, Republic of Korea. 
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FIGURES 

Fig. 1. Decay Pj^ distributions of a) 531 KINK (muonic one prong) events, and b) 276 VEE 
(muonic two prong) events. The points with error bars are the experimental data, and the histograms 
are the Monte Carlo simulations. The vertical dotted lines show the decay Pj^ limits of the main 
semi-leptonic modes. In (a), the solid, dotted and dashed histograms represent the final fit result, 
the sum of the backgrounds, and the — > l^v^ signal, respectively. In (b), the solid line represents 
the inclusive spectrum of muons from D° decay. 
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TABLES 



Table 1. Summary of offline selection criterion and cut requirements. 



Offline selection 
Flight length cut 
Decay Pj^ cut 

Vertex position cut (only for VEE) 



Pr^Beam > 0.8 GeV/c 

FL > 500 lira 

P Tll > 0.28 GeV/c 

Decay vertex in emulsion bulk 



Table 2. Cross sections used for Monte Carlo background calculations. 



a(D ± ; x F > 0) 8.66 ± 0.46 ± 1.96 fib/nucleon [ 8 ] 

a(D°; x F > 0) 22.05 ± 1.37 ± 4.82 /ib/nucleon [ 8 ] 

<r(Df; x F > 0) 5.1 ± 1.3 ± 1.3 /ib/nucleon [ 9 , 18 ] 

a(Af; x F > 0) 15.3 ± 1.7 ± 2.6 /ib/nucleon [ 10 ] 



Table 3. Summary of published measurement of pseudoscalar decay constant /d s . 



Experiment /d s (MeV) 

E653 (This work) 194 ±35 ±20 ±14 

WA75 f 1 ] (1993) 238 ± 47 ± 21 ± 43 

CLEO [ 2 ] (1994) 344 ± 37 ± 52 ± 42 

BES [ 3 ] (1995) 430±igg ± 40 
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